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HH , Abstract A model of low-frequency quasi-periodic oscillations (LFQPOs) of black hole 

r~| ■ X-ray binaries (BHXBs) is proposed based on the perturbed magnetohydrodynamic (MHD) 

Qh. equations of accretion disk. It turns out that the LFQPOs frequencies of some BHXBs can be 

fitted by the frequencies of the toroidal Alfven wave oscillation corresponding to the max- 
\Zi \ imal radiation flux. In addition, the positive correlation of the LFQPO frequencies with the 

j/j ■ radiation flux from accretion disk is well interpreted. 
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00 ■ 1 INTRODUCTION 

o. 

(^ , It is well known that the quasi-periodic oscillations (QPOs) have been widely observed in some X-ray bi- 

^J ■ naries, whose compact object is either a neutron star or a black-hole, and QPOs play an essential role as 

^—v \ a potentially important tool for studying the strong gravitational field and understanding the physical pro- 

^\i ■ cesses of X-ray states (Done, Gierlinski & Kubota 2007). Among ^20 black hole X-ray binaries (BHXBs), 

the low-frequency QPOs (LFQPOs) with frequency ranging from a few mHz to 30Hz have been detected 
on one or more occasions for 14 systems in the hard state and the steep power-law (SPL) state (McClintock 
& Remillard 120061 1. Almost all the LFQPOs have some common properties, remaining relatively stable 
^^ , features for days or weeks usually. 

^ ' A number of theoretical models have been proposed to explain the physical mechanisms of LFQPO. 

Tagger and Pellat ( 1999^ suggested that LFQPOs observed in low-mass BHXBs can be interpreted by 
accretion-ejection instability. Titarchuk & Osherovich (120001 ) considered that the LFQPOs are caused by 
the global disk oscillation in the direction normal to the disk, and these oscillations arise from the grav- 
itational interaction between the central compact object and the disk. O'Neill et al. (1 20111 1 proposed that 
LFQPOs are related to the quasi-periodic behavior in global MHD dynamos. Cabanac et al. ( 120101 ) pre- 
sented that an oscillating corona due to magneto-acoustic wave propagating in the corona, and produce 
multiple QPOs. Kato (2008) proposed that the one-armed c-mode low frequency oscillations of disk are 
one of possible candidates of LFQPOs based on a resonantly-excited disk-oscillation model. In addition, 
relativistic precession model was firstly presented by Stella and Vietri (1998), they suggested that LFQPOs 
are the result of some modulation of Lense-Thirring precession. Later, Schnittman et al. (2006) and Ingram 
& Done et al. (12009120 12] l developed this model, and successfully explained the variability properties of 
BHXBs. However, there has been no consensus on the physical nature of LFQPOs. 

On the other hand, it has been detected that LFQPO frequencies exhibit a strong positive correlation 
with disk flux from observations. For example, LFQPO frequencies of XTE J1550-564 (33 observations 
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during 1998.9-1999.4), GRO J1655-40 (38 observations during 2005.2.17- 2005.3.6) and H1743-322 (20 
observations during 2003.3.28-2003.5.22) display a roughly linear relation with disk flux, when they are in 
the hard and intermediate states with frequencies in the range l-7Hz, 0.01-2.5Hz and 2-5Hz, respectively 
(Remiflard & McCHntock ,2006. McClintock et al. |2009l Sobczak et al. l2000l Shaposhnikov et al. .2(J07] l. 
Similar quasi-linear relations of LFQPO frequencies with disk flux have been found for GRS 1915H-105 (32 
observations during 1996. 10.29-1998. 10.7) in the range 1- 15Hz (Markwardt, Swank & Taam 1999"; Sobczak 
et al. 2000; Muno, Remillard & Morgan 200JJ and forXTE J1748-288 ( 92 obs ervations during 1998.7.13- 
1998.9.26) in the range 20-30Hz (Revnivtsev, Trudolyubov & Borozdin 120001 1. From the duration time of 
the above outbursts, we can find that this positive correlation between the LFQPO frequencies and disk 
flux is on long timescales. It is noted that the QPO observed in XTE J1550-564 during its 1998 outburst 
shows a correlation between absolute rms amplitude and mean source flux over timescales shorter than 
^3ks (Heil,Vaughan & Uttley 120111) . so this relation also holds on short timescales. However, the origin of 
the correlation between LFQPO frequency and disk flux of BHXBs remains elusive. 

If magnetic field is taken into account in accretion disk, a torsional Alfven wave can be generated by 
the rotational dragged of space (Koide et al. 120021) . We consider a toroidal magnetic field existing in a 
rotational accretion disk, in which Alfven wave oscillation propagates along toroidal magnetic field lines 
due to the perturbation of radial velocity of the accreting matter Thus the Alfven wave oscillation will 
influence the transformation of angular momentum and the radiation flux from the inner disk. Although 
some researchers suggested that the QPOs of low mass BHXBs arise from Alfven wave (Zhang 120041 Shi 
& Li l2010t Shi l2011l ). the association of Alfven wave oscillation with LFQPOs has not been discussed. 

In this paper, we adopt a binary system consisting of a Kerr black hole surrounded by weakly magne- 
tized relativistic thin disk to derive the Alfven wave frequency based on the perturbed MHD equations, and 
the LFQPO frequency is fitted by the frequency of Alfven wave propagating in a specific circular orbit of 
the accretion disk, which corresponds to the maximum radiation flux. In addition, we propose a thin disk 
model with some corona floating above the disk to well interpret the positive correlation of the LFQPO 
frequencies with disk flux. 

This paper is organized as follows. In section 2, we present a description of our model, and derive 
the formula for Alfven wave frequency and radiation flux. In section 3, we fit the LFQPO frequencies of 
several BHXBs, and fit the positive correlation of LFQPO frequency with the disk flux for two BHXBs, 
XTE J1550-564 and GRO 1655-40. Finally, in section 4, we discuss the results obtained in our model. 

2 MODEL DESCRIPTION 

2.1 Alfven Wave Oscillations in Accretion Disk 

We consider a geometrically thin, optically thick, non-self-gravitating perfect fluid disk, which is magne- 
tized and isothermal, and rotating around a BH. It is assumed that the magnetic pressure is much less than 
the total pressure(p„iag ^ p), which is a good approximation to a thin Keplerian disk (Li 120021) . The disk 
dynamics is governed by the ideal MHD equations given as follows. 

|^+V.(pv)=0, (1) 

^ + v- Vv+-Vp + V$+-; — Bx(VxB)==0. (2) 

ot p Anp 

The magnetic field B satisfies the induction equation in the MHD approximation, and it reads 

-^ - V X (v X B) = 0, (3) 

and 

V • B = 0. (4) 

The quantities p, v and p in equations (1) and (2) denote mass density, velocity and pressure of plasma, 
respectively. The Pseudo-Kerr potential (Mukhopadhyay |2002| l $ is adopted to simulate the effects of 
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general relativity, and it reads 

V., = - cV^2a.VF + «;)^ (5) 

where r is the distance between the plasma to center of BH (r = R/Rg), and Rg = GM/c^ is the 
gravitational radius. The quantities M,G and c denote respectively the BH mass, gravitational constant and 
speed of light, and J and a* = J /{GM/c^) represent the BH angular momentum and the dimensionless 
spin, respectively. 

The Keplerian angular velocity in the frame of Pseudo-Kerr potential is derived by Shafee et al. (I2008I I 
as follow, 

c3(r^-2a,Vf + a^) 

'' GMr2[v^(r-2) + a*]- ^' 

The perturbed physical quantities in MHD equations can be written as 

p = Po + p', P = Po + p', V = vo + v', B = Bo + B', (7) 

where the subscript '0' and superscript '" denote the equilibrium and perturbation values, respectively. 
Assuming that all the perturbations are small enough, i.e., po ^ p' , Po ^ p'^ "^o ^ ^'^ Bq ^ B', we can 
neglect the products of them in second and higher orders. Substituting equations (4) and (7) into equations 
(l)-(3), we have the linearized MHD perturbation equations as follows, 

do' 

^+Po(Vv') + (Vp')-vo) = 0, (8) 

^+vo-(Vv') + ^Vp' + -^x(VxB')=0, (9) 

dt po 47rpo 

QTD/ 

-- - V X (V X Bo) = 0, (10) 

ot 

From the vertical equilibrium assumption of accretion disk, the half-height(i/) of the disk can be written as 
(Mukhopadhyay 2003) . 

H = c,R'/^F-^/^^cjnk, (11) 



where Cg — \/p/p is the sound speed. 

Incorporating equations (8)-(10), we have the equation for the velocity perturbation as follows, 

^^ + vo • V( — ) + c2v(V • v') + VA • V X [V X (v' X va)] = 0, (12) 

where va = ,^° is Alfven velocity, being defined as 

VA = ^^ - v^c, = ^/^H^k, (13) 

where /? = Pmag/p is defined as the ratio of magnetic pressure to total pressure. 

In order to fit LFQPO frequencies by invoking Alfven wave propagating in an accretion disk, we express 
the perturbed velocity as follows, 

v' = v'e'^*'-?-"*), (14) 

where ui, ^ and k are the perturbation frequency, the displacement vector and the wavenumber vector, re- 
spectively. In cylindrical coordinates the wavenumber vector is k = krV + k^(p + fc^z, where kr, k^, kz 
represent the radial, azimuthal and vertical components, respectively. Substituting equation (14) into equa- 
tion (12), we have 



,2,/ 



- wfcvo • v' + {cl + 4)(k • v')k + (va • k) • [(va • k)v' - (va • v')k - (k • v')va] = 0. (15) 
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Considering that Alfven wave is a transverse wave, transporting along the magnetic field line, we have 

k||Bo||vo J- v', and equation (14) is simplified as 

k^v\ - tj2 = 0, (16) 

where the symbols || and _L denote parallel and vertical directions, respectively. Thus the angular velocity 
of the perturbation is related to the Alfven velocity va as follows, 

W — kvA = ky^Cs- (17) 

For a thin magnetized accretion disk, Armitage and Natarajan (1 19991 1 presented that the dominant field 
component is toroidal with saturation occurring when Pmag ^ P- So we assume that the unperturbed 
flow is axisymmetric, and only toroidal magnetic field exists in equilibrium state in the accretion disk,i.e.. 
Bo — B^ip. Thus a small perturbation of radial velocity of fluid due to accretion of black hole gives 
rise to the oscillation of the magnetic field B^p, and results in the Alfven wave oscillation transporting 
in the toroidal direction in the disk. For the characteristic wavelength X ^ R, the toroidal wave number 
k^ ~ 2tt/R (Shi. l2010] l. and the angular velocity of the perturbation can be written as 

UJ = k^^Cs^27T—^nk. (18) 

H 

Inspecting equation (18) we find that u is much less than the Keplerian angular velocity ilk in the thin 
disk due to the disk scale height H being much less than disk radius R, and this perturbation frequency 
provides a possibility for fitting LFQPO frequencies of BHXBs. 

2.2 Relation between LFQPO Frequency and Perturbation Frequency 

First of all, we intend to clarify the relation between LFQPO frequency and the perturbation frequency. 
The disk angular velocity can be regarded as the Keplerian angular velocity, i.e., il = ilk, provided that 
the radial magnetic force can be neglected. This result can be realized if only toroidal magnetic field exists 
without vertical electric current in the thin disk. 

The rate of energy generation per unit area of one side of the disk is given by (Shakura & Sunyav 119731 
Novikov & Thorne \fW3[ Gierliriski et al. [T999l Shafee et al. l2008] l. 

Q{R) = -RHap— = ^^ — , (19) 

where M and a are respectively the mass accretion rate and viscosity parameter, ilj^ and flms denote 
respectively the angular velocity at R and inner edge of disk. The radius R„is of the innermost stable 
circular orbit (ISCO) of accretion disk varies with a* in Pseudo-Kerr potential, Mukhopadhyay (120021) 
presented that r^s (rms — Rms/Rg) satisfies the foUowing equation 

- 3at + 14ai^^Aw + (r™s - 6)r^, + 6a^rli^{rras + 2) - 2alr„^sir„,s + 11) = 0. (20) 

It is easy to find that ISCO moves to the BH with the increasing a* as shown in Fig. 1. Thus the local 
radiation flux F{r) can be written as: 

where we have r = R/Rg, and 

g{r,r„is,a^) = (2r^V^W - r^V^r^- r^a* + 2a*y^^- a^\/^ + rf„^V^- 2rf, 



ms V 



+rl^sa* - 2a^^/r3r,ns + alV^){12V^a^ - 3r* + IGa^r - Wa^V^ - 7r^al 
+2r'' - 4al^)/[^{V^ - 2V^ + a,fi^/^, - 2^1— + a,)]. (22) 
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Fig. 1 The curves of local radiation flux F{r) versus r with the given values of a, 
lOAf0,M = 0.1MEdd,MEdd is the Eddington accretion rate.) 



(M = 



By using equation (20) and (21), we plot the curves of local flux F{r) versus r as shown in Fig.l, and 
we find that the radiation flux is dominantly produced in the inner disk, and it varies non-monotonically 
with r, attaining its peak value at r^ close to ISCO. In addition, as shown in Fig.l, we find that the peak 
value of F{r) is greater, and the location of r^ is closer to ISCO for a greater BH spin. 

Inspecting equations (18)-(21), we can find that the radiation flux is directly caused by the angular mo- 
mentum transport. For magnetized thin accretion disk, numerical simulations have shown that the angular 
momentum transport is dominated by Maxwell stresses (Hawley et al. 119951 Brandenburg et al. 119951 Stone 
et al. 1 19961 1, and the magnetic contribution to viscosity parameter a exceeds fluid stresses by an order of 
magnitude, so that a ~ amagnetic (Armitage and Nataraian |1999l l. So in low mass X-ray binaries the high 
energy X-ray radiation is generated mainly from the interaction of the plasma with the magnetic field, and 
the change of the magnetic field could modulate the X-ray flux (Shi 2011). Thereby the oscillations of the 
toroidal Alfven wave may lead to the QPOs. Since the strongest influence of the perturbation on luminosity 
corresponds to the maximum radiation flux F^ax {r = rd), the perturbation frequency at r^ can be regarded 
as the LFQPO frequency. 

3 FITTING LFQPOS OF BHXBS 

In this section we intend to fit the LFQPOs of several BHXBs based on the above relation between LFQPO 
frequency and the perturbation frequency. In addition, the positive correlation between LFQPO frequency 
and disk flux can be well interpreted based on our model. 

We consider a geometrically thin, optically thick accretion disk with corona floating above given by 
Gierliriski et al. (I1999I I. In the inner region of weak magnetic fields, the gas pressure and magnetic pressure 
can be neglected, and the total pressure is 



P — Prad 



4cr 
3c 



Tr 



(23) 



where a is Stefan-Boltzmann constant, and the quantity Tc is the temperature of the central disk, which 
is determined by the equilibrium between the radiation cooling in the vertical direction and the energy 
generated by the viscous dissipation. We assume that a fraction q of the total energy generated by the 
viscous process is emitted from the disk, and the remainder is dissipated in the corona. This energy equation 
is written as (Frank, King & Raine .2002i Gierlihski et al. 1 19991 1 



4cr 

37 



T^ = qF{r) 



(24) 
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Table 1 Fitting LFQPO Frequencies of X-ray BHBs. 

source i^qpo rn^BH "» P = Pmag/p m 

GROJ1665-40 0.1-28 6.2 0.7 0.00086-0.24 

XTEJ1550-564 0.1-10 9.1 0.34 0.01 0.0023-0.23 

GRS 1915+105 0.001-10 15.0 0.975 0.00007-0.07 

4U 1543-47 7 9.4 0.8 0.073 

Notes: "Remillard & McClintock ( |2006l >; ''Narayan & McClintock ( |20TB 

The parameter r in equation (24) is the opacity for electron scattering, and it reads 

r = pHar/mp, (25) 

where ctt and nip denote the Thomson cross-section and proton mass, respectively. Incorporating equation 
(11) with equations (23)-(25), we obtain 



mpcirtkY 



H - ,„ JA .^ F{t). (26) 



VQPo = :^ == 5-1 X 10''—- 3 ^^ '^^ ^ j^ 7^^5(^rf,''ms,a,). (29) 



The Eddington accretion rate MsdcL can be written as 

MEdd = 2TTRjnsmpc/r]aT, (27) 

where r/ is the accretion efficiency, whose value is 0.1 for black hole. Substituting equations (6), (21) and 
(27) into equation (26), we have 

H = 3.7 X lO^mBHmq—- — JT^— 51^: "Tms, «*), (28) 

[Zasfy'r — r — a^) 

where we define itibh = M/Mq and rn = AI/MEdd- 

Combining equations (6), (18) and (28), we have the perturbation frequency as a function of disk radius 
r, and we derive the expression for LFQPO frequency by setting r = r^ as follows, 

-^ == 5 1 X io5 ^"^g^""(V^~^V^+"*) , 
27r ' rriBH ^d(2a*V^ ~ ^d ~ "*) 

From equation (29) we find that the LFQPO frequency is proportional to accretion rate m, but it is 
inversely proportional to the BH mass m. Based on the above results we can fit some LFQPO frequencies 
of several BHXBs with the given BH mass by adjusting the mass accretion rate and the intensity of magnetic 
field /3, and the main results are summarized as follows. 

First, the theoretical values of the LFQPO frequencies are in accordance with the observed ones with 
appropriate values of m and /3. We assume that the dissipation in corona is neglected (q = 1), and /3 = 
Pmag/p = 0.01, and the fitting results are listed in Table 1. 

Second, the observed strong positive correlation between the LFQPO frequencies and the disk flux can 
be interpreted very well based on our model. 

We assume that the disk locally emits a blackbody spectrum, of which a fraction 1 — psc is scattered in 
the corona, and the disk luminosity is (Gierlihski et al. 119991 ) 

L, = PscqvMc^ = — —, (30) 

cos 2 

where i and D are the inclination of the disk and the source distance to the observer, respectively. The 
quantity (1 — psc) is the fraction of the disk emission which is not scattered by the corona, and Fg is 
the observed disk flux. Combining equations (29) and (30), we obtain the relation between the LFQPO 
frequency and the observed disk flux Fg as follows. 

VQPO = 2.5 X 10-33^i^i^ : 7,7 ^^ ^2 ^g{rd, r^s, a* • (31) 
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Table 2 LFQPO Frequencies Fitting Parameters of XTE J1550-564 and GRO J1655-40. 



Input 

niBH 



D{kpc) 



Ps 



Output 



rd 



Fitting parameter x^ (d-of) 

P = Pmag/P 



XTEJ1550-564 9.1 0.34 4.38 75° 0.15 4.83 7.64 0.0105 ± 0.0003 44.9(32) 
GROJ1665-40 6.3 0.7 3.2 70° 0.1 3.39 5.33 0.0159 ± 0.0005 61.2(37) 

Notes:The values of input parameters im fin.a*.D.i) are adopted from Naravan & McClintock ( 1201 It . 
The quantities of r,„s and ra are computed by using equations (20) and (21). 




0.2 0.4 



Fig. 2 Fitting the relation between the disk flux and LFQPO frequency for XTE J 1550-564. 
The observation data are taken from Fig. 2 Of SobcZak et al. (2000), while the flux errors are 
corrected based on McClintock et al. (2009). 



From equation (31) we find that LFQPO frequency has a positive correlation witii flie observed disk 
flux. For XTE J1550-564 (outburst in 1998) and GRO J1655-40 (outburst in early 2005). The values of 
Psc are about 15% (Sobczak et al. 2000) and 10% (Shaposhnikov et al. ,2007) . respectively, and they keep 
constant in hard state. Therefore their disk fluxes have a linear relation with LFQPO frequencies. 

For a given source, the values of rms and rd can be obtained by resolving equations (20) and (21), and 
the values of the parameters, such as itibh, a^,,D, i, and psc)are available in the literatures. Substituting all 
the parameters into the equation (31), we can derive the LFQPO frequency corresponding to disk flux Fs 
and the magnetic field parameter f3 as follows. 



l^QPO = fW)Fs, 



(32) 



where f{f3) is the function of /3. The measurement values of LFQPO frequency and disk flux are Vi and 
Fi (i — 1,2,- ■■ ,N, the number of the observation data), of which the errors are e^i and ej-i.We adopt the 
Nukers' estimate (Tremaine et al. 12002) based on the following minimizing 



TV 



2 _ 






fWFi. 



PWh 



(33) 



Fi 



Using the Nukers' estimate, we can derive the best fitting parameter (3 for XTE J1550-564 and GRO 
J1655-40, which are presented in Table 2, and the values of reduced x^ per degree of freedom are less than 2 
for both sources, indicating that LFQPO frequency has approximate linear relation with disk flux as shown 
in Figs.2 and 3 for XTE J1550-564 and GRO J1655-40, respectively. 
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Fig. 3 Fitting the relation between the disk flux and LFQPO frequency for GRO J1655-40. 
{iTiBH = 6.2, j3 ~ 0.6, D = 3.2kpc, i = 70°. The observation data are taken from Fig. 15 of 
Shaposhnikov et al. (2007). 



4 DISCUSSION 

In this paper, we propose that the LFQPOs in BHXBs can be interpreted by invoking the toroidal Alfven 
wave oscillations located at the disk radius with the maximal radiation flux. It turns out that the LFQPO 
frequencies of several BHXBs can be well fitted based on our model. In addition, the positive correlation 
between the LFQPO frequencies and the disk flux of XTE J1550-564 and GRO J1655-40 are also wefl 
fitted. According to the argument given in sections 2 and 3, we find that the positive correlation between 
LFQPO frequency and the disk flux can be understood, because both the perturbation frequency of the 
toroidal Alfven wave and LFQPO frequency increase with the increasing accretion rate. It is noted that this 
correlation only holds in hard stats, when the disk fraction keeps constant. Thus our model provides an 
explanation for the fact that LFQPOs are primarily observed in the power law part of spectra of BHXBs 
(Zycki & Sobolewska|2005]l. 

In this paper, we adopt a thin disk with some corona floating above, which is different from the corona 
interior to a truncated disk as given by Done, Gierlinski & Kubota (2007). The truncated disk is not com- 
patible with our model because it is difficult to determine the location of maximal radiation flux. If Alfven 
wave of truncated disk is used to interpret QPOs, a certain appropriate location of Alfven wave propagating 
must be found out, e.g. Shi ( 1201 11 1 used Alfven wave oscillation at transition radius to explain HFQPOs of 
low mass X-ray binaries. 

Although LFQPOs are interpreted successfully by invoking perturbed Alfven wave oscillation, it seems 
difficult to explain high Frequency QPOs (HFQPOs) of BHXBs. This could imply that the physical origins 
of LFQPOs and HFQPOs are different. HFQPOs of BHXBs have been successfully explained by the mag- 
netic hot spots model (Wang, et al. 2003 2005 ) and resonance mode (Abramowicz and Kluzniak 20011 
Huang, et al. 1201 Oi l, and the HFQPO frequencies are sensitive to BH mass and spin, while they seems not 
directly related to the accretion rate and disk flux. It is noted that HFQPOs may be triggered by the insta- 
bility of accretion disk oscillation (Tagger et al. 2006, 2001). We also notice that the instability occurs in 
the toroidal Alfven wave oscillation, and the oscillation frequency is comparable to the HFQPO frequency, 
provided that the accretion rate and the magnetic field are great enough. This result motivates us to explore 
the relation of HFQPOs and the instability of Alfven wave oscillation in our future work. 
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